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Electron-electron interactions generally reduce the low temperature resistivity due to the screen- 
ing of the impurity potential by the electron gas. In the weak-coupling limit, the magnitude of this 
screening effect is determined by the thermodynamic compressibility which is proportional to the 
inverse screening length. We show that when strong correlations are present, although the com- 
pressibility is reduced, the screening effect is nevertheless strongly enhanced. This phenomenon is 
traced to the same non-perturbative Kondo-like processes that lead to strong mass enhancements, 
but which are absent in weak coupling approaches. We predict metallicity to be strongly stabilized 
in an intermediate regime where the interactions and the disorder are of comparable magnitude. 

PACS numbers: 72.15.Qm, 71.55.Jv 



Transport in disordered metals has been studied for 
many years, and substantial theoretical and experimental 
understanding has been achieved 0] in the case of weak 
disorder and in the regime of weak electron-electron inter- 
action. Much less is known about situations with strong 
electronic correlations; here, most research has concen- 
trated on clean systems. 

Recent experiments on two-dimensional (2D) electron 
gases in zero magnetic field have led to considerable 
renewed interest in electronic systems close to metal- 
insulator transitions (MITs). In these systems, well- 
defined metallic behavior (positive temperature coeffi- 
cient of resistivity, dp/dT > 0) has been observed in 
the low density regime, and is characterized by a sur- 
prisingly large (up to a factor of ten) drop of resistivity 
at low temperatures. Here, the electron-electron inter- 
actions represent the largest energy scale in the problem 
0; this is emphasized by recent reports of substantial 
mass enhancement from several complementary experi- 
ments Q. 

Some simple microscopic mechanisms that can produce 
such a resistivity drop relate to temperature-dependent 
screening 0, 0, El Q of the random potential. This 
effect obtains even in the simplest Hartree-Fock (HF) 
treatment of interactions, which represents the basis for 
the standard (Lindhard) screening theory, and applies 
equally well to both short-range ja] and long-range Q 
forces. In this picture, the screening length is inversely 
proportional to the thermodynamic compressibility of the 
system, which therefore controls the magnitude of disor- 
der screening. This observation immediately brings into 
question the relevance of such mechanisms in the regime 
of strong correlation, in particular close to interaction- 
driven MITs. Here, one expects substantial mass en- 
hancements, but at the same time a significant decrease 
of the compressibility ||. Then, if applicable, standard 



screening theory (as in HF) would predict weak disorder 
screening precisely where mass is enhanced - in contrast 
to experiments, where the resistivity drop persists in that 
region (see Fig. 3, below). 

In this letter, we examine the screening of the impu- 
rity potential by focusing on a model where a reliable 
and controlled treatment of strong correlations is avail- 
able. This is possible within dynamical mean-field the- 
ory (DMFT) which is formally exact in the limit of 
large coordination. To investigate the regime of strong 
correlation, we examine a disordered Hubbard model in 
the vicinity of the Mott transition, a model which has 
recently been argued 0, 0, 0, 0] to provide an appro- 
priate description of interaction effects near the 2D-MIT. 
Our results demonstrate that: (1) If DMFT equations are 
solved within the HF approach, we reproduce the results 
of standard screening theory in qualitative and even semi- 
quantitative detail; in this picture, screening is strongly 
suppressed close to the Mott transition. (2) A more ac- 
curate solution agrees with HF results far from the Mott 
transition, but finds diametrically opposite results in the 
regime of strong correlation. Here, while compressibility 
is reduced, both the disorder screening and the effective 
mass are strongly enhanced. (3) The enhanced screening 
strongly stabilizes metallic behavior in the intermediate 
regime where the disorder and interactions are of com- 
parable magnitude. 

We consider a disordered Hubbard model described by 
the Hamiltonian 



(1) 



Here tij are the hopping matrix elements, c and c 1 are 
fermionic creation and annihilation operators, n = c'c is 
the number operator, and a labels the spin projection. U 
represents the Hubbard on-site repulsion, and the disor- 
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der is introduced by random site energies £j, as specified 
by a distribution function P(si). 

Disorder renormalization (screening). Within DMFT 
, a quasiparticle is characterized by a self-energy func- 
tion S(w), which is assumed to be purely local (momen- 
tum independent). In a random system 0], this quan- 
tity though still local, is now site-dependent = 
Si) which explicitly depends only on the correspond- 
ing local site energy Ei The renormalized disorder 
potential (as seen by the quasiparticle at the Fermi en- 
ergy) can thus be defined by 



Vi{£i) 



+ (uj = 0) — (I — 6(1, 
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where (i is the chemical potential. These renormalized 
energies are defined with respect to a reference energy 
6(i chosen such that vl = J dei P(ei)vi(ei) = 0, i.e. that 
their site-average vanishes. In a case of particle- hole sym- 
metry 8(i — 0. Since there are no vertex corrections 
within DMFT, the T — dc resistivity d E3 depends 
only on the variance of this renormalized disorder, viz. 
p ~ vf. The self-energies E^w) must be calculated by 
solving an ensemble of Anderson impurity problems sup- 
plemented by a self-consistency condition 0, 0] . 

Weak coupling (Hartree-Fock) solution. In the HF ap- 
proximation for the Hubbard model, £i(w) = U (nj) ; 
here (• ■ • ) represents the quantum average (for a given 
disorder configuration). Within DMFT, the local occu- 
pation (rii) depends only on the local (renormalized) site 
energy, and for moderate disorder Q we find 



\ = [i + u Xii ]- 2 ., 
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(3) 



where Xti — ~d jdsi is the local compressibility. [lq 
This result, although based on a local approximation, 
proves to provide qualitative and even semiquantitative 
agreement with more standard screening theory. Just as 
in Refs. [J0-II1- f° r reasonable values of the interaction 
strength (U < bandwidth), the screening of the random 
potential cannot be very large (resistivity drop by a factor 
of two at mostQ). 

Strongly correlated regime. We approach aT = Mott 
transition at half-filling. In general, the DMFT equa- 
tions cannot be solved in closed form, but extensive work 
H has shown that most qualitative and even quantita- 
tive features of the T — solution can be reproduced 
using simple analytical approximations. Here, we use 
the four-boson mean-field method of Kotliar and Ruck- 
enstein (KR),Q which is equivalent to the well-known 
Gutzwiller variational approximation, but can be readily 
generalized to disordered systems. This approach pro- 
vides a parameterization of the low-energy (quasiparti- 
cle) part of the local Green function which, in our case, 
takes the form (h = c = 1) 

GiM = ~ ~ Z \ a t r . (4) 
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FIG. 1: DMFT phase diagram of the Hubbard model with 
random site energies. The Mott insulator can be suppressed 
by sufficiently strong (bare) disorder W > U. Also shown is 
an estimate of the regime where Anderson localization effects 
are important, as obtained by comparing the Fermi energy 
Ef to bare disorder W (dotted line) or screened disorder W 
(dashed line). Localization is strongly suppressed by correla- 
tion effects in the intermediate regime where the disorder is 
comparable to the on-site repulsion U. 



Here, the local quasiparticle weight Zi = 
2 [I - {el - d?) 2 ] 1 (e, + dtf [1 - (ef + d?)] , as well 
as parameters e^, di and el are site-dependent 
quantities, determined by the KR equations 0] 



dZi i 



■q-^o X! &-( ljJ n )G l (uj n ) = Zi (fi + gj - Si) e t 



de, (3 
dZi I 



,•;,/ p X! A(uj n )G l (uj n ) = Zi (U-fj,-Si +£i) d{, 
I^G i K) = ^ i (l-e 2 + d 2 ). 
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Half- filling can be enforced by the requirement e 2 = df. 
Finally, a self-consistency condition 12| determines the 
"hybridization function" A(uj) = A (u> — Yi av (u))) de- 
scribing the environment of a given site. Here A (ui) = 
u> + fi— [G (u>)]~ , G (u>) being the lattice Green's func- 
tion corresponding to no disorder and U = 0. The "av- 
erage" self-energy E „(a>) is defined via the disorder- 
averaged Green's function G(uj) = J dei P(si)Gi(u)): 
G(lu) = G {ijj — S otJ (a;)) . Thus S aw (w) = uj + (i — A(ui) — 

[G(w)] . It follows that 5(i = J dei P(ei)si/Zi, and 
Vi = e l jZ i - Sfi. 

Phase diagram. These equations can be easily solved 
for an arbitrary band structure and distribution of disor- 
der, but such details do not affect the qualitative form of 
the solution, which proves to depend only on the presence 
or absence of particle-hole symmetry. We first examine 
the particle-hole symmetric situation, and as an illus- 
tration we concentrate on a model with a semi-circular 



band of width At at half-filling, and site energies uni- 
formly distributed in the interval (— W/2, W/2). The re- 
sulting T = phase diagram, as obtained from a full 
numerical solution, is shown in Fig. 1. The Mott insu- 
lating phase is completely suppressed for W > U, since 
the disorder tends to fill in the Hubbard-Mott gap. The 
phase boundary separating the correlated metal and the 
Mott insulator (full line) is identified by the simultane- 
ous vanishing of the quasiparticle weights Zi on all lattice 
sites. All quantities display simple critical behavior close 
to this phase boundary, the form of which can be ana- 
lytically obtained for weak disorder, but which proves to 
remain qualitatively correct along the entire critical line. 
To second order in £j, we find Zi — Zq(\ + Cef), where 
Z = Zi(si = 0) = 2(1 — U/Uc) linearly goes to zero at 
the critical interaction U C (W) = U°(l+ fCef). This cor- 
responds to the well-known effective mass enhancement 
(m* ~ 1/Z) in the strongly correlated regime. Here the 
constant C = 10/(£/ c ) 2 , and U° = 64£/3tt for the consid- 
ered band structure. 

Scattering rate. To determine the effects of correla- 
tions on transport, we calculate the scattering rate 1/r = 
-2ImE a „(cj = 0) = -21mG o (0)vf. We find that scat- 
tering is strongly reduced near the Mott transition, cor- 
responding to correlation-enhanced screening. To show 
this analytically, we note that in the critical regime the 
spectral weight corresponding to A(w) is of the order 
of Z Q . We expand Eq. J7J) with respect to e,-/^: To 
leading order, ii/Zi ~ — e 2 + d 2 ~ O(Zi). Therefore 
EijZi = Vi ~ Zi —> as U — > U C (W), and we conclude 
that random site energies are perfectly screened at the 
metallic side of the Mott transition |l6j . The same con- 
clusion remains valid even for an asymmetric distribution 
of disorder, in which case particle-hole symmetry is re- 
stored as the transition is approached. For weak disorder 
we obtained a simple formula for the scattering rate close 
to the Mott transition 

r--t[W c \ Uc)\ V () 

Typical numerical results obtained for disorder ranging 
from weak to strong are shown in Fig. 2(a). We find that 
disorder is strongly screened even relatively far from the 
transition, for U/U C (W) > 0.5. 

Particle-hole asymmetry. For a particle-hole asymmet- 
ric model (generally appropriate for realistic materials), 
we find that for U —> U C (W) the energies Vi remain fi- 
nite, although still much smaller than W. The disor- 
der screening nevertheless remains very strong, which 
persists even in the case of fairly strong particle-hole 
asymmetry and strong disorder. Typical numerical re- 
sults are shown in Fig. 2(b), where the full line corre- 
sponds to the symmetric semicircular density of states, 
and the dashed line is obtained for the particle-hole asym- 
metric lattice described by the Green function G (ui) = 




u/u c 



FIG. 2: Scattering rate normalized with the noninteract- 
ing value, (a) From the lower to the upper curves: W = 
0.1, 1,2,4. (b) Results close to U c for the particle-hole sym- 
metric lattice (full line), and the asymmetric lattice (dashed 
line), with W = 1. The inset shows the density of states in 
these two cases. 

\u> — t 2 G {uj + a)] . In this plot, the asymmetry pa- 
rameter a is chosen to be 0.1 5D (strongly asymmetric 
lattice), where D is the bandwidth. For weak disor- 
der and moderate particle-hole asymmetry we were able 
to obtain a quantitative estimate of the screening ef- 
fect, and we find that close to the transition the scat- 
tering rate approaches a very small asymptotic value 
1/r ~ (5)i) 2 ef ~ W 4 , instead of 1/r ~ W 2 as for weak 
interactions. In addition, for any realistic lattice at half 
filling, (5/x is a small number (< 0.1), explaining the small- 
ness of the scattering rate. 

Breakdown of conventional theory. In the strongly- 
correlated regime our DMFT results are in sharp discrep- 
ancy with results obtained within Hartrcc-Fock, as shown 
in Fig. 3. For small U, both methods give similar results, 
but closer to the transition HF theory (dashed line) pre- 
dicts a reduced disorder screening, while full DMFT (full 
line) shows that the screening remains strongly enhanced. 
The reduction of screening found in HF reflects the de- 
crease of the compressibility near the Mott transition. 
This is a result of a Stoner instability in the magnetic 
HF solution, which sets in for U — l/po(0) = irD/2 as 
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FIG. 3: Normalized scattering rate as a function of U, from 
the full DMFT solution (full line), and the corresponding 
Hartree-Fock (HF) approximation (dashed line). For moder- 
ate interaction both methods predict the same screening, but 
diametrically opposite results are obtained in the strongly- 
correlated regime, where DMFT predicts enhanced screening, 
while a strong suppression is obtained within HF theory. 

a precursor to a gap opening at the transition [ItJ. We 
emphasize that the direct relation between compressibil- 
ity and screening as in Eq. © is a general feature of 
weak-coupling approaches, and its applicability is seri- 
ously limited in the strongly-correlated regime. 

Enhanced screening as "Kondo pinning". In the 
DMFT approach || that we use, the solution of the full 
Hubbard model is mapped to solving an ensemble [Lj] of 
auxiliary Kondo- Anderson impurity problems. Accord- 
ingly, the approach to the Mott transition can be de- 
scribed as the decrease of the local Kondo temperature, 
corresponding to the reduction of the local quasiparti- 
cle weight. In this language, the renormalized energy 
level Vi can be identified as the position of the Kondo 
resonance, which is well-known to "pin" to the Fermi 
energy in the Kondo limit Zi — > 0. We can thus in- 
terpret the surprising enhancement of disorder screen- 
ing in the strongly-correlated regime as reflecting the 
non-perturbative Kondo physics captured by our DMFT 
method, but not by standard weak-coupling theories. 
This mechanism is very closely related to the Kondo- 
enhancement of resonant tunneling through quantum 
dots, as observed in recent experiments 0]. O ur discus- 
sion makes it clear why site randomness is strongly sup- 
pressed, but also indicates that if additional hopping ran- 
domness |l2l is introduced, the same mechanism would 
not apply, as we have verified by explicit calculations. 
In realistic systems we expect the disorder to be strongly 
but not perfectly screened even in the vicinity of the Mott 
transition. 

Screening and localization. The DMFT approach is 
too simple to describe Anderson localization effects, 
which cannot be neglected for strong enough random- 
ness. Nevertheless, it is interesting to estimate the disor- 



der strength necessary for localization. In the absence of 
interactions, localization is expected to set in when the 
disorder scale W is comparable to the kinetic (i.e. Fermi) 
energy, as indicated by a dotted line in Fig. 1. However, 
we have shown that correlations lead to strong screening, 
with a renormalized disorder scale W ~ (vf) 1 ' 2 <C W, 
which we can numerically compute for any U and W, and 
analytically in several limits. In particular, in the atomic 
limit (E F 0), we find W ~ (1 - U/W) 3/2 . In the 
presence of interactions, thejDnset of localization should 
be estimated by comparing W to Ep, and the resulting 
boundary is shown by a dashed line in Fig. 1. Inter- 
estingly, the metallic phase is found to be strongly sta- 
bilized by screening in the intermediate regime W ~ U. 
Of course, such interplay of correlation and localization 
should be studied in more detail by extensions of DMFT 
which can explicitly incorporate the localization effects 
[l9) . but this fascinating issue remains a challenge for 
future work. 
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